Compressed entanglement testing 
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We introduce a sequence of numerical tests that can determine the entanglement or separability 
of a state even when there is not enough information to completely determine its density matrix. 
Given partial information about the state in the form of linear constraints on the density matrix, 
the sequence of tests can prove that either all states satisfying the constraints are entangled, or 
that there is at least one separable state that satisfies them. The algorithm works even if the 
values of the constraints are only known to fall in a certain range. If the states are entangled, 
an entanglement witness is constructed and lower bounds on entanglement measures and related 
quantities are provided; if a separable state satisfies the constraints, a separable decomposition is 
provided to certify this fact. 



Entanglement is one of the central features in quan- 
tum information processing (QIP). It has been identified 
as a key ingredient in many useful QIP tasks such as tele- 
portation, quantum key distribution, superdense coding 
and quantum computation [l[ . It is also remarkable that 
even when a complete description of a quantum state is 
given (in the form of a density matrix), it can be ex- 
tremely difficult to computationally decide whether such 
state is entangled or not. This is due to the fact that 
this problem (usually called "the separability problem" ) 
is known to be NP-Hard 0. The problem is even more 
difficult when we consider experimental tests of entangle- 
ment, since measurements may not provide a full descrip- 
tion of the state, and when they do (such as in quantum 
state tomography [3||) the reconstructed density matrix 
may be unphysical (i.e., not positive semidefinite (PSD)). 

A key problem with important practical applications is 
to determine the entanglement characteristics of a state 
when only a limited amount of information is available. 
If this information comes from measuring a set of ob- 
servables, it takes the form of a set of linear constraints 
on the elements of the density matrix. This set of con- 
straints may also come from other more theoretical con- 
siderations. In this letter we will introduce a sequence 
of numerical tests that can decide whether all states that 
satisfy a set of linear constraints are entangled, or if there 
is at least one separable state that satisfy those same con- 
straints. If the states are shown to be entangled, the al- 
gorithm constructs an entanglement witness that certifies 
this fact for all such states and provides lower bounds on 
certain entanglement measures and related quantities. If 
a separable state satisfying the constraints exists, the al- 
gorithm finds it and provides a separable decomposition 
as a proof. 

Consider a situation in which we are given partial in- 
formation of the state of a quantum system in the form of 



L linear constraints on the elements of its density matrix 



Tr[pM z ] 



I = 1. 



(1) 



where the operators Mi are arbitrary. Our goal is to de- 
termine if all the states satisfying these constraints are 
entangled, or if there is a separable state that satisfies 
them. The constraints in ([1]) arc nothing but a linear 
system of equations for the elements of the density ma- 
trix p. If this system is incompatible it means that these 
constraints do not describe a physical state. If the sys- 
tem is invertible, then the Hcrmitian matrix p can be 
completely determined from the equations, and once we 
have an explicit expression we can check if it corresponds 
to a state (i.e., it is PSD and normalized), and then ap- 
ply any of a battery of separability criteria. In particu- 
lar, if we apply the combined PPT symmetric extension 
(PPTSE) criterion introduced in |4( and its dual intro- 
duced in Q , we know that we can conclusively decide the 
entanglement properties of the state (although it may be 
computationally very expensive for some states). But the 
situation that is the most interesting (and typically more 
common) corresponds to the case in which the linear sys- 
tem ([T|) is undcrdctermined, and we do not have enough 
information to uniquely define the state. This situation 
corresponds naturally to being able to measure only the 
expectation value of a limited number of observables (the 
operators Mi in ([T]) are then Hermitian matrices). We 
will show that in this case, the linear system defines an 
affine subspace in the state of Hermitian matrices, and 
the PPTSE criterion and its dual can be applied to either 
prove that all states in that affine subspace are entangled, 
or to show that a separable state exists that satisfies ([1} . 

Let us start with the linear system ([1]). The most gen- 
eral solution p of this system can be written as 



Dj, 



part 



(2) 
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where (P art is a particular solution of ([T]) (i.e., 
Tr[pP art Mi] = mi, I = 1,...,L), the matrices {^ a )} 
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form a basis of the subspace of solutions of the homoge- 
neous system (i.e., Tr[^ a )M ; ] =0, I = 1, . . . , L), D K is 
the dimension of this subspace, and y a are real variables. 
Note that p part is just a Hermitian matrix and not nec- 
essarily a state (i.e., it need not be PSD or normalized). 
The question then reduces to whether there are values 
of the real variables y a such that the resulting Hermi- 
tian matrix is a normalized, separable state. If there are 
not, then all the states of the form that is all nor- 
malized PSD Hermitian matrices satisfying (JIJ must be 
entangled. 

In order to test the separability of © we will use the 
PPTSE criterion Q. We define a PPT symmetric ex- 
tension of p to k copies of subsystem A as a state p 
in Hb such that, (i) p = Tr^-i [p], (ii) p is 

symmetric under exchanges of copies of subsystem A, 
and (Hi) p has positive partial transposes for any bi- 
partite arrangement of the subsystems A and B. Since 
any separable state in Ha ® Hb can be written as 
p = ^2 PiV^ij {^i I ® \4>i) {4>i 1 1 it has such an extension given 
by P = ilPiWi){i>i\® k ® |0»)v&|- For each value of k, 
the non-existence of a PPTSE provides a sufficient (but 
not necessary) condition for entanglement. In the limit 
k — y oo the condition becomes necessary. The practi- 
cal value of this approach is that searching for such ex- 
tensions or proving their impossibility can be cast as a 
semidefinite program. 

A semidefinite program (SDP) is a type of convex opti- 
mization problem that has a broad range of applications 
and has been widely applied in quantum information. An 
SDP has both a primal and a dual form. A typical SDP 
in its primal form reads 



minimize c T x 
subject to F + J2i x i F i h 0, 



(3) 



where c is a given vector, x = (x\, ... ,x n ), and Fq and 
Ft are some fixed Hermitian matrices. The inequality in 
the second line means that the affine combination of the 
F matrices must be positive semidefinite. The minimiza- 
tion is performed over the vector x, whose components 
are the variables of the problem. The dual of this SDP 
takes the form 



maximize — Tr[i*oZ] 
subject to Z >z 0, 

Tr[FiZ] = a, 



(4) 



where the dual variables are the components of the ma- 
trix Z . 

To test the entanglement of a state of the form ©, we 
can apply the PPTSE criterion for any value of A:. We will 
present in detail how this works for k — 2 (the general 
case is straightforward) . So we need to check if, for some 
values of the variables y a , the resulting matrix is PSD, 

normalized and has a PPTSE. Let {of}^, {af }fg 1 be 
bases for the spaces of Hermitian matrices that oper- 
ate on Ha and Hb, of dimensions <1a and ds respec- 
tively, such that they satisfy Tr^cr^er^] = ax$ij and 



Tr[crf-] = 8n (where X stands for A or B), and ax 
is some constant. Then we can expand p in the ba- 
sis {af (g> of}, and write p = V / ; pijaf (g> of , with 
t^a^Tilpaf ® of]. In the same way, we can 



Pij = 

expand the extension p in H® 2 <g> Hb as 



^PikjWf 



i<k 



^2 P kk 3 a k 
kj 



(5) 



where we made explicit use of the swapping symmetry 
between the two copies of A. To satisfy the condition 
that p is an extension of p, we need to impose Tr^ [p] = p. 
This implies puj — py. From ^ we have that = 

Pij Tt + J2a=i VaPij^ i which fixes some of the components 
of the extension ([&]) . We then have 

P = £ of + of ® at ® of } + 



i>l 



D K I 

■ y* J2 ® a f +a i® °i ® of }+ 

a=l \ « 
\i>\ 

- 51 p*feiK A ® ® of + CT fe ® ffA ® °f ) + 

i>fe>2 

■^PfcfciOfc ®<^®of. (6) 



? ft 

fe>2 



If we define a vector of variables x = (y, pikj ) (with 
2 < k < i < d 2 A . 1 < j < d? B ), we can see that the 
most general form of the extension ([6]) has the form 
Go + '^2 i XiGi, where the expressions for Go and Gi can 
be easily extracted from it. 

The first condition we need to impose on this exten- 
sion is that it represents a state, i.e., that it is PSD 
and normalized. The normalization condition can al- 
ways be assumed to be contained in the set of linear 
equations |T|), by adding another constraint with M = 1 
and expectation value equal to 1. Requiring that the 
extension is PSD means imposing the linear matrix in- 
equality (LMI) G + X)i x * G i h 0. And finally, im- 
posing the positivity of the partial transposes requires 
two more LMIs, namely G A + J2i x iGj A h and 
G o B + E l %i G J B h (due to the swapping symme- 
try, these are the only two independent partial trans- 
poses). We can combine these three LMIs into a single 
one by defining matrices F = Go © Gj. A © Gq b and 
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Fi = d © Gf A ®Gj B (a block diagonal matrix is PSD if 
and only if all of its blocks are PSD). So searching for a 
PPTSE of a state of the form © corresponds to a SDP of 
the form ([3J with c = (0, . . . , 0). If there are values of y a 
such that |2]) is separable, then there must exist values of 
p ik j, (2 < k < i < d\, 1 < j < d%) such that the SDP is 
feasible (because separable states always have PPTSE). 
But if the SDP is infeasible, it means that there is no set 
of values y a for which the resulting state p has a PPTSE, 
and hence all states of the form © must be entangled. 

We can take this approach a little bit further and con- 
sider the case in which the expectation values of the op- 
erators Mi (the RHS of ([I])) are known only approxi- 
mately. Assume that instead of ((T|) we have Tr [pMj] S 
[mj" m ,mj" ax ], I = 1, ...,L. Now consider the set of 
matrices {r p : Tr [r p M;] = 5 p i}. We can use these matri- 
ces to write any particular solution of the linear system 
Tr[pAf/] = zi as pP art = J2f=i z i T i- Combining this with 
@ we have 

L D K 

P = J2^n + Y,ya^ a \ (7) 

1=1 o=l 

as the most general solution of Tr[pMj] <E 
[m,y iin ,m 1 j nax ], I = 1,...,L, provided that 

2; e [m| Tlin ,mJ 7loa! ]. If we apply the PPTSE crite- 
rion to ([7]) as before, we will obtain a linear combination 
of matrices representing the PPT symmetric extension 
of a state satisfying this set of equations. We just 
need to once again construct the required LMIs to 
impose the positive semidcfinitcncss of the extension 
and its partial transposes, and solve the resulting SDP 
satisfying the constraints zi € [mj™ n , mj" ax ]. These 
constraints can be imposed by another LMI, namely 
diag(zi - m^ m , m^ ax -zi,...,z L - m^ m , m% ax -z L )> 
0, showing that constraining the range of the variables 
does not change the SDP structure. 

Another useful feature of the PPTSE criterion is that, 
if the primal SDP is infeasible (i.e., there is no separable 
state satisfying the constraints), the dual SDP provides 
a certificate of this fact in the form of an entanglement 
witness [f|. Let us recall that an entanglement witness 
is a Hermitian operator whose expectation value is posi- 
tive on all separable states but negative on some entan- 
gled states. It provides a proof of the entanglement of a 
given state. Entanglement witnesses are a particular case 
of the separating hyperplane theorem (or Hahn-Banach 
theorem) of convex geometry: if two closed convex sets 
are disjoint and one of them is compact, there is a hyper- 
plane that separates them (i.e., the sets are in opposite 
sides of the hyperplane). In the context of checking sep- 
arability of linearly constrained states, the convex sets 
in question are the set of separable states and the affinc 
subspace spanned by all the solutions of the linear sys- 
tem {!]). If these two sets are disjoint, it means that no 
separable state satisfies ([T]); on the other hand, the sep- 
arating hyperplane theorem assures us that there is an 
entanglement witness that can certify the entanglement 



of every state of the form Q . 

A slight modification of the SDP allows us to com- 
pute bounds on certain useful entanglement quantities. 
Let us add an extra variable t and an extra term to the 
LMI: F + J2i x i F i + tl >: 0. This problem is equiva- 
lent to the original if we minimize t: if the original is 
feasible, t opt < 0, and t opt > if it is infeasible. But 
then d 2 A d,Bt pt is a lower bound on the minimum amount 
of the maximally mixed state we need to add to a state 
satisfying fl} to make it separable (in 2 ® 2 and 2 ® 3 
this bound is tight). This is known as the random ro- 
bustness of entanglement j.7], R r (p), and quantifies how 
robust the entanglement is against white noise. It also 
provides a lower bound on a geometric measure of quan- 
tum discord |8|. 

The entanglement witness constructed from the dual 
SDP can also be used to quantify the entanglement of the 
states satisfying |T]). Any entanglement measure that can 

be expressed as E(p) = max{0, ~^j^Tr[Wp]} with M 
a subset of entanglement witnesses, is referred to as wit- 
nessed entanglement Q. The set A4 determines which 
particular measure this expression represents. Several 
well-know measures are of this form, such as the best sep- 
arable approximation BSA(p), the negativity Af(p), and 
the concurrence C(p). Clearly, any W € Ai that satisfies 
Tr[Wp] < provides a lower bound to E(p). In partic- 
ular, the quantities E n ^ m {p) — max{0, ^g™ 11 Tr[VKp]} 
(n, m > 0) with M n . m = {W : -nl ^ W < ml} a sub- 
set of entanglement witnesses, are entanglement mono- 
tones, and satisfy E n<m (p) — > nBSA(p) when m oo, 
where BSA(p) is the best separable approximation to 
p [9(. Since E n ^ m (p) is obviously monotonically increas- 
ing with m (for fixed n) and any entanglement witness Z 
must be in some A4 n ,rm Tr[Zp] provides a lower bound on 
BSA(p), which is an entanglement measure. This anal- 
ysis is just an illustration of the connection between the 
PPTSE criterion and entanglement measures, and does 
not pretend to give the best bounds possible. 

Each test in the PPTSE hierarchy provides a sufficient 
but not necessary condition for entanglement. The hi- 
erarchy is complete in the limit: any entangled state 
is guaranteed to be detected by one of the tests. But 
a separable state will pass all tests, leading to a non- 
terminating algorithm. Fortunately, a dual approach was 
developed by Navascues et al.,@, that applies a sequence 
of tests that can certify separability in a finite number of 
steps (although that number can be very high for some 
states). Geometrically, the PPTSE hierarchy of tests 
works by monotonically approximating the cone of sep- 
arable states from the outside with a sequence of cones 
associated with states having PPT symmetric extensions 
to a certain number of copies of one of the subsystems. 
The dual approach in Q constructs a similar approxima- 
tion to the cone of separable states, but from the inside: 
it provides sufficient (but not necessary) conditions for 
separability. By interleaving the two sequences of tests 
we can, in a finite number of steps, determine if a state is 
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entangled (and give an entanglement witness as a proof), 
or separable (and provide an explicit separable decompo- 
sition). We will now briefly describe the test in [5j and 
show how it can be used to determine if there is a state 
of the form ([2]) that is separable. 

Let Sp be the set of states in Ha Ki'Hb that have a 
PPT symmetric extension to N copies of A. In [l(| it 
was shown that a small perturbation in Hb makes these 
states separable. More precisely we have that Sp = {(1 — 

£n)uab + £nwa® I s : uab € Sp } satisfies S^ c S, for 
all N, where S is the set of separable states in Ha <8> Hb 
, wa = Tt b [ujab], and e N = d B /(2(d B - l))min{l - 
x : Pjfr^+i m ° d2) (x) = 0}, with Pt 0) (x) the Jacobi 
polynomials. Since Sp — > Sp for N going to infinity, and 
S C S]? for all N, we have that S^ -> S (N -> oo). This 
result can be easily transformed into a SDP that tests if a 
given state is separable [5] . The same approach applies to 
states defined by (0) , and so we can construct a sequence 
of tests that is guaranteed to find a separable state that 
satisfies (TTJ) if one exists. Moreover, an explicit separable 
decomposition will be provided in that case (the actual 
form can be found in [10|). 

Let us use a simple example to illustrate the 
power of this approach. Consider a system that 
produces two photons and we want to determine if 
they are entangled in the polarization basis. One 
possible approach is to do quantum state tomogra- 
phy. This can be accomplished by measuring the 
16 observables given by [1_J] fii ® fij(i,j = 0,1,2,3) 
with Ao = \H){H\ + \V)(V\, Ai = \H)(H\,fa = 
\D){D\,ji 3 = \R)(R\, with \D) = (\H) - \V))/V2 
and \R) = (\H) - i\V))/V2. Note that these op- 
erators are all positive on pure product states and 
so they are good candidates to be entanglement 
witnesses. Applying our test we find an exam- 
ple in which only four of these observables are 
enough to prove entanglement: if Tr[(Ai <8> Ai)/°] G 
[0.48, 0.5], Tr[(Ai ® £_)p] G [0.24, 0.25], Tr[(A 2 ® fa)p] € 
[0.48,0.5], Tr[(p 3 ®p, 3 )p\ G [0,0.02], (note that all expec- 
tation values are non negative) , there is an entanglement 
witness Z = 0A343\HH}(HH\ + 0.3977\HV)(HV\ + 
0.234(\VH)(VH\ + \VV)(VV\) + {(0.0658 + 
iO. 1583) (\HH}(VH\ + \HV)(VV\ + \VH)(VV\) + 



h.c.} + {~0.2242\HH)(VV\ + 0.0925\HV)(VH\ + h.c}, 
such that Tr[Zp] < —0.0168 for all states satisfying the 
linear constraints above. Moreover, Z G Aii.i, so this 
result provides a lower bound on the best separable 
approximation, i.e., BSA(p) > 1.68 x 10~ 2 . The primal 
SDP also computes a lower bound on the random 
robustness, R r (p) > 8 x 0.0168 = 0.1344. Additional 
information about the state can improve these bounds. 
For example, if we add Tr[(/ti ® p. 3 )p] G [0.24, 0.25] to the 
constraints we now obtain a new entanglement witness 
Z' such that Tr[Z'p] < -0.021, which translates to 
BSA(p) > 2.1 x 10~ 2 , and R r (p) > 0.168 (the MATLAB 
code used is available online from the author [12(). 

An interesting point to consider is when does this com- 
bined approach fail, i.e., it does not produce a definite 
answer after a reasonable number of steps. If we look at 
the PPTSE part of the method (proving entanglement), 
we can see that if the affine space defined by © intersects 
sets Sp with N large, even if the states are actually en- 
tangled it will take at least N steps to prove this fact. On 
the other hand, if we consider the dual approach and (J2|) 
intersects the set of separable states but does not inter- 
sect Sp for N small, the procedure will take a long time 
to provide the required separable decomposition. Since 
and Sp approach S from the outside and the inside 
respectively, the most challenging situation for this pro- 
cedure occurs when ([2]) is "almost" tangent to the set 
of separable states, either intersecting it or not. Fail- 
ure of the approach to provide an answer in a reasonable 
number of steps suggests that stronger constraints on the 
state may be needed. 

In summary, the algorithm described in this letter al- 
lows us to determine the entanglement character of states 
even when only partial information about them is avail- 
able. It also provides bounds on measures of entangle- 
ment and other related quantities. This tool can thus 
be very useful to experimentally determine if a state is 
entangled when the number of observables that can be 
measured is limited. 

I would like to thank Sergio Boixo, Daniel Lidar and 
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tions that helped to refine the approach presented in this 
work. 
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